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Abstract 



Scattering of electromagnetic (EM) waves by many small particles 
(bodies), embedded in a thin layer, is studied. Physical properties 
of the particles are described by their boundary impedances. The 
thin layer of depth of the order 0{a) with many embedded in it small 
particles of characteristic size a, is described by a boundary condition 
on the surface of the layer. The limiting interface boundary condition 
is obtained for the effective EM field in the limiting medium, in the 
limit a — ^ 0, where the number M{a) of the particles tends to infinity 
at a suitable rate. 
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1 Introduction 

It is known (see, e.g., [2], [9]) that the light propagation through diffraction 
gratings may exibit strong resonances at certain frequencies. This is useful in 
applications. In this paper we study electromagnetic (EM) wave scattering 
by many small impedance particles Dm, 1 < m < M, M = M(o), embedded 
in a thin layer of the depth h{a) ~ a, where a is the characteristic dimension 
of a small particle. The shape of the particles may be fairly arbitrary, not 
necessarily spherical. 
We assume that 





ka < 1, 
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where k is the wavenumber, and d{a) is the distance between neighboring 
particles. 

The thin layer is located on a smooth surface S. The permettivity eo, 
conductivity o"o > 0, e' = eo + ^ and permeability /xq of the space are 
assumed constants, k'^ = uj'^e'fio, to is the frequency. 

For example, one may assume that S is the plane X3 = 0, but our 
arguments are valid for an arbitrary smooth S. The M particles on S are 
distributed according to the following law: in any open subset A of 5" there 
are AA(A) particles, where 

M{^) = ^ f N{s)ds[l + o{l)], a^O, (2) 

N(s) > is a continuous function, vanishing outside of a finite domain 
O C 5 in which small particles (bodies) Dm are distributed, k G (0, 1) is a 
number, and the boundary impedances of the small particles are defined by 
the formula 

_ h[sm) c n ("i^ 

where Sm G 5* is a point inside m— th particle Dm-, Re > 0, and /i(s) is 
a continuous function vanishing outside 0. 
We can choose k and h{s) as we wish. 

Denote by [E, H] = E x H the cross product of two vectors, and by 
{E, H) = E ■ H the dot product of two vectors. 

The impedance boundary condition on the surface Sm of the m— th par- 
ticle Dm is = (^m[H^ , N], where E^ (H^) is the tangential component of 
E (H) on Sm, and N is the unit normal to Sm, pointing out of Dm- We 
define -E* = [N, [E, N]] = E — N{E, N). This corresponds to the geometrical 
meaning of the tangential component of E, and differs from the definition 
E^ = [N, E] that is used sometimes. 

In this paper we use the methodology, developed in |2j and some results 
from [7]. 

2 EM wave scattering by many small particles 

Electromagnetic (EM) wave scattering problem consists of finding vectors 
E and H satisfying the Maxwell equations: 

V X E = iu;fioH, V x H = -icoeoE in D :=R^ \\J^^^Dm, (4) 
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the impedance boundary conditions: 



[N, [E, N]] = Cm[H, N] on l<m<M, (5) 

and the radiation conditions: 

E = Eo + VE, H = Ho + VH, (6) 

where ve and vh satify the radiation condition, (m is the impedance, defined 
in (l3|), is the unit normal to Sm pointing out of D^, Eq,Hq are the 
incident fields satisfying equations ([T|) in all of . One often assumes that 
the incident wave is a plane wave, i.e., Eq = f e*'^"'^, £^ is a constant vector, 
a G S*^ is a unit vector, S'^ is the unit sphere in M^, a ■ £ = 0, ve and 
Vh satisfy the radiation condition: — ikv) = oil) as r := |x| — )• oo, 

k = u^eofio- 

By impedance we assume in this paper a constant, Re Cm > 0, or a 
matrix function 2x2 acting on the tangential to 5m vector fields, such that 

Re(CmE\E')>0 ^E'eTm, (7) 

where is the set of all tangential to Sm continuous vector fields such 
that Divi?* = 0, where Div is the surface divergence, and E*' is the tan- 
gential component of E. Smallness of means that ka <^ 1, where 
a = 0.5 maxi<m,<M diamDm- 

Lemma 1. Problem ([4]) -([7]) has at most one solution. 

Lemma [His proved in j7]. 
Let us note that problem (Il])-([7]) is equivalent to the problems ([8]), ([9]), ([6]), 
d?]), where 

V xV X E = k'^E in D, H = ^ ^ ^ , (8) 

[N, [E, N]] = -^[V X E, N] on Sm, 1 < m < M. (9) 

This is the impedance boundary condition (see, e.g., [Ij, p. 301.) The 
expression [A^, [E, N]] = E — {E, N)N is the tangential component of the 
field E on the surface Sm, is the unit normal to Sm pointing out of Dm- 

Thus, we have reduced our problem to finding one vector E[x). If E{x) 
is found, then H = 

Let us look for E of the form 

I^I ft %]i^x y I 

E = Eo+y^Vx g{x,t)am{t)dt, g{x,y) = ^ (10) 

^-^ q 47r X — y\ 

m=l I 



4 



where t G Sm, dt is an element of the area of S'm, and (Tm{t) G T^. This E 
for any continuous crm{t) solves equation ^ in D because Eq solves ([8|) and 

VxVxVx / g{x,t)a.m{t)dt = VV - V X / g{x,t)am{t)dt 

- V^V X / g{x,t)am{t)dt 

J Sm 

= fc^V X / g{x,t)am{t)dt, x G D. 

^ Sm 

(11) 

Here we have used the known identity divcurlE = 0, valid for any smooth 
vector field E, and the known formula 

-V^g{x,y) = k^g{x,y) + d{x-y). (12) 

The integral Jg ^((x, t)(Tm(t)di satisfies the radiation condition. Thus, for- 
mula pOj) gives E{x) that solves problem (HI), ([9]), and satisfies the radiation 
condition, if (Jm{t) are chosen so that boundary conditions ([9]) are satisfied. 

Define the effective field Ee{x) = E^{x) = E^\x,a), acting on the 
m— th body Dm'- 

Eeix):=E{x)-V X f gix,t)am{t)dt:=Ei"'\x), (13) 

where we assume that x is in a neigborhood of Sm- However, Ee{x) is defined 
for ah X G M^. 

Away from S, the field Ee{x, a) tends to a limit E(x) = Ee{x) as a — )• 0, 
and Ee{x) is a twice continuously differentiable function away from S, see 
[7]. To derive an integral equation for am = CTmit), substitute 

E = Ee + Vx I g{x,t)am{t)dt 

J Sm 

into boundary condition Q, use the known formula (see, e.g., [3]) 
[iV,Vx / g{x,t)am{t)dt]^ = j [Ns,[V sg{x,t),am{t)]]dt ±^^, (14) 

J Sm Sm 

where the -(+) signs denote the limiting values of the left-hand side of (|14|) 
as X — )■ s from D {Dm)-, and get 

(Tm{t) = AmfJm + fm, 1 < m < M. (15) 
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Here Am is a linear Fredholm-type integral operator, and fm is a continu- 
ously differ entiable function. 

Let us specify Am and fm- One has 

/^ = 2[iV„/e(s)], Ms):=[Ns,[E,{s),Ns]]-^[V xEe,Ns]. (16) 
Condition ^ and formula (fH|) yield 

/e(s) + ^[amis),Ns] + [ / [N„ [Vsgis,t),amit)]]dt,N,] 



Srr 



(17) 



[VxVx / g{x,t)am{t)dt,Ns]\^^s = 

Using the known formula V x V = graddiv — V^, the relation 
VxVx- g{x,t)am{t)dt = Vx {-Vtg{x,t),amit))dt 

= Vx j g{x,t)BWam{t)dt = 0, 
where Div is the surface divergence, and a consequence of formula ([12 



(18) 



-Vt g{x,t)am{t)dt = e g{x,t)am{t)dt, x ^ S, (19) 

one gets from (fT7|) the following equation 

[N,,am{s)] + 2fe{s) + 2Bam = 0. (20) 

Here 



Bam-=[ [Ns,[Vsg{s,t),amit)]]dt,Ns]+Cmiu}eo[ g{s,t)amit)dt, N,]. 

Sm S m 

(21) 

Take cross product of Ns with the left-hand side of ()20p and use the formulas 
Ns ■ (Tmis) = 0, fm ■= fm{s) := 2[Ns, fe{s)], and 

[Ns,[Ns,am{s)]] = -am{s), (22) 

to get from (pO|) equation (fTSj) : 

a„(s) = 2[Ns, fe{s)] + 2[iV„ := A^a^ + (23) 



where AmCTm = '^[Ns,Bam\- The operator Am is hnear and compact in the 
space C{Sm), so that equation (j23|) is of Predholm type. Therefore, equation 
([23]) is solvable for any fm G Tm if the homogeneous version of (p3|) has only 
the trivial solution am = 0. In this case the solution am to equation (j23|) is 
of the order of the right-hand side fm, that is, 0{a~'^) as a — )• 0, see formula 
()16p . Moreover, it follows from equation (j23p that the main term of the 
asymptotics of — )• does not depend on s G Sm- 

Lemma 2. Assume that am £ Tm, cfm £ C{Sm), end am{s) = Amf^m- Then 
am = 0. 

Lemma [2] is proved in [7]. 
Let us write (fTOl) as 



M M „ 

S(x) = ^o(a;)+^[V^5f(x,x„),Qm]+^ Vx / {g{x,t)-g{x, Xm))crmit)dt, 

m=l m=l •'•^"i 

(24) 

where 

Qm-= [ Om{t)dt. (25) 



Since am = 0(a ''), one has Qm = 0(a^ '^). We want to prove that the 
second sum in (j24p is negligible compared with the first sum. One has 



[V,5(x,x^),Q™]| < O (^max|^,^|^ ©(a^^'^), (26) 

aO ^max - 



j2:=|Vx / {g{x,t)-g{x,Xm))crm{t)dt\<aO[uia.-}^{^,^\]0{a^ 



d^' d 



(27) 

and 

J2 

ji 

Thus, one may neglect the second sum in ([26l) . and write 



O (^max|^,fca}) ^ 0, ^ = ^(^)' a ^ 0. (28) 



M 



E{x) = Eoix) + ["^x9ix, Xm),Qm] (29) 



m=l 



with an error that tends to zero as a — t- 0. Let us estimate Qm asymptoti- 
cally, as a — )■ 0. Integrate equation (p3|) over 5^ to get 



Qm = 2 I [Ns, fe{s)]ds + 2 f [Ns,Bam]ds. 



(30) 
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It follows from (fT6l) that 



[Ns, fe] = [Ns, E,] - -^[Ns, [V X Ee, Ns]]. (31) 

If Ee tends to a finite limit as a — )■ 0, then formula (I3ip implies 

[Ns, fe] = OiU) = O (^^) , a ^ 0. (32) 

By Lemma[2]the operator (/ — Am)~^ is bounded, so am = O (^), and 

Qm = (a^-'^) , a ^ 0, (33) 



because integration over Sm adds factor O(a^). As a — )• 0, the sum (p9]) 
converges to the integral (see Lemma 1) 

E = Eo + Vx / g{x,s)N{s)Q{s)ds, (34) 
Js 

where N{s) is the function from ([2]), and Q{s) is the function such that 

Qm = Qixm)a^-\ (35) 
The function Q{y) can be expressed in terms of E: 

Q{s) = -^h{s)iVxE)is), (36) 

see [7]. Here the factor ^ appears if are balls. Otherwise a tensorial 
factor Cm, depending on the shape of Sm, should be used in place of The 
factor Cm is defined by the formula /^^ V x Ee{s)ds = a^'CmS x E(,{xm), 
where Xm G S is a point in Dm- 

Thus, equation takes the form 

E{x) = Eo{x) ^Vx / g{x,s)h{s)N{s)V X E{s)ds. (37) 

3w/io Js 

It follows from equation (j37p that the limiting field £'(x) satisfies equa- 
tion dHI) away from 5, and a transmission boundary condition on S: 

[N,,E.{s) - E+{s)] = -^h{s)N{s)V X E{s). (38) 
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3 Conclusions 

It is proved that a distribution of many small impedance particles in a thin 
layer on a smooth surface S can be described by a transmission boundary 
condition ()38p . This condition shows that the equivalent surface currents 
on S are calculated analytically in the limit o — )• in terms of boundary 
impedance function h and the distribution density function N{s). Therefore, 
these currents can be controlled. 
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